A vertex‐centered linearity‐preserving discretization of diffusion problems on polygonal meshes

This paper introduces a vertex‐centered linearity‐preserving finite volume scheme for the heterogeneous anisotropic diffusion equations on general polygonal meshes. The unknowns of this scheme are purely the values at the mesh vertices, and no auxiliary unknowns are utilized. The scheme is locally conservative with respect to the dual mesh, captures exactly the linear solutions, leads to a symmetric positive definite matrix, and yields a nine‐point stencil on structured quadrilateral meshes. The coercivity of the scheme is rigorously analyzed on arbitrary mesh size under some weak geometry assumptions. Also, the relation with the finite volume element method is discussed. Finally, some numerical tests show the optimal convergence rates for the discrete solution and flux on various mesh types and for various diffusion tensors. Copyright © 2015 John Wiley & Sons, Ltd.     

A FULLY DISCRETE IMPLICIT-EXPLICIT FINITE ELEMENT METHOD FOR SOLVING THE FITZHUGH-NAGUMO MODEL

This work develops a fully discrete implicit-explicit finite element scheme for a parabolic-ordinary system with a nonlinear reaction term which is known as the FitzHugh-Nagumo model from physiology. The first-order backward Euler discretization for the time derivative, and an implicit-explicit discretization for the nonlinear reaction term are employed for the model, with a simple linearization technique used to make the process of solving equations more efficient. The stability and convergence of the fully discrete implicit-explicit finite element method are proved, which shows that the FitzHugh-Nagumo model is accurately solved and the trajectory of potential transmission is obtained. The numerical results are also reported to verify the convergence results and the stability of the proposed method.     

Asymptotic error distribution for the Euler scheme with locally Lipschitz coefficients

In traditional works on numerical schemes for solving stochastic differential equations (SDEs), the globally Lipschitz assumption is often assumed to ensure different types of convergence. In practice, this is often too strong a condition. Brownian motion driven SDEs used in applications sometimes have coefficients which are only Lipschitz on compact sets, but the paths of the SDE solutions can be arbitrarily large. In this paper, we prove convergence in probability and a weak convergence result under a less restrictive assumption, that is, locally Lipschitz and with no finite time explosion. We prove if a numerical scheme converges in probability uniformly on any compact time set (UCP) with a certain rate under a global Lipschitz condition, then the UCP with the same rate holds when a globally Lipschitz condition is replaced with a locally Lipschitz plus no finite explosion condition. For the Euler scheme, weak convergence of the error process is also established. The main contribution of this paper is the proof of $\sqrt{n}$ weak convergence of the normalized error process and the limit process is also provided. We further study the boundedness of the second moments of the weak limit process and its running supremum under both global Lipschitz and locally Lipschitz conditions.     

Convergence issues in using high‐resolution schemes and lower–upper symmetric Gauss–Seidel method for steady shock‐induced combustion problems

This paper reports numerical convergence study for simulations of steady shock‐induced combustion problems with high‐resolution shock‐capturing schemes. Five typical schemes are used: the Roe flux‐based monotone upstream‐centered scheme for conservation laws (MUSCL) and weighted essentially non‐oscillatory (WENO) schemes, the Lax–Friedrichs splitting‐based non‐oscillatory no‐free parameter dissipative (NND) and WENO schemes, and the Harten–Yee upwind total variation diminishing (TVD) scheme. These schemes are implemented with the finite volume discretization on structured quadrilateral meshes in dimension‐by‐dimension way and the lower–upper symmetric Gauss–Seidel (LU–SGS) relaxation method for solving the axisymmetric multispecies reactive Navier–Stokes equations. Comparison of iterative convergence between different schemes has been made using supersonic combustion flows around a spherical projectile with Mach numbers M = 3.55 and 6.46 and a ram accelerator with M = 6.7. These test cases were regarded as steady combustion problems in literature. Calculations on gradually refined meshes show that the second‐order NND, MUSCL, and TVD schemes can converge well to steady states from coarse through fine meshes for M = 3.55 case in which shock and combustion fronts are separate, whereas the (nominally) fifth‐order WENO schemes can only converge to some residual level. More interestingly, the numerical results show that all the schemes do not converge to steady‐state solutions for M = 6.46 in the spherical projectile and M = 6.7 in the ram accelerator cases on fine meshes although they all converge on coarser meshes or on fine meshes without chemical reactions. The result is based on the particular preconditioner of LU–SGS scheme. Possible reasons for the nonconvergence in reactive flow simulation are discussed.Copyright © 2012 John Wiley & Sons, Ltd.     

A finite volume scheme preserving extremum principle for convection–diffusion equations on polygonal meshes

We propose a nonlinear finite volume scheme for convection–diffusion equation on polygonal meshes and prove that the discrete solution of the scheme satisfies the discrete extremum principle. The approximation of diffusive flux is based on an adaptive approach of choosing stencil in the construction of discrete normal flux, and the approximation of convection flux is based on the second‐order upwind method with proper slope limiter. Our scheme is locally conservative and has only cell‐centered unknowns. Numerical results show that our scheme can preserve discrete extremum principle and has almost second‐order accuracy. Copyright © 2017 John Wiley & Sons, Ltd.     

The parameter uniform numerical method for singularly perturbed parabolic reaction–diffusion problems on equidistributed grids

This article presents the study of singularly perturbed parabolic reaction–diffusion problems with boundary layers. To solve these problems, we use a modified backward Euler finite difference scheme on layer adapted nonuniform meshes at each time level. The nonuniform meshes are obtained by equidistribution of a positive monitor function, which involves the second-order spatial derivative of the singular component of the solution. The equidistributing monitor function at each time level allows us to use this technique to non-linear parabolic problems. The truncation error and the stability analysis are obtained. Parameter–uniform error estimates are derived for the numerical solution. To support the theoretical results, numerical experiments are carried out.     

Air Flow Computation around an Automated Guided Vehicle

An upwind finite element scheme for the incompressible viscous flow at a high Reynolds number was proposed by the fourth and fifth authors. The scheme has the potential to approximate the advection term in third-order accuracy. We apply it to a two-dimensional non-stationary analysis of airflows around an Automated Guided Vehicle (AGV), which starts with constant acceleration, runs at a constant speed and stops with constant deceleration. The results are at least qualitatively good and compatible with experimental ones.